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1. Introduction 



In Monte Carlo simulations of unquenched lattice QCD, a positive-definite fermion determi- 
nant is required to ensure the convergence of the Markov chain. For Hybrid Monte Carlo (HMC) 
simulation [jl[, it also requires that the pseudofermion action is positive-definite such that the pseud- 
ofermion field can be simulated with a Gaussian noise. For QCD with two degenerate flavors, this 
is guaranteed provided that det(D t ) = det(D). The challenge of simulating one fermion flavor with 
HMC is that it is nonuivial to write down a positive-definite pseudofermion action with determi- 
nant exactly equal to det(D). So far, there are several strategies to simulate one fermion flavor with 
HMC, e.g., to approximate D 1 by a polynomial of even degree or to approximate (D^D) -1 / 2 
by a rational polynomial function In this paper, we construct a positive-definite pseudofermion 
action with determinant exactly equal to det(D), for Wilson and domain- wall fermions respectively. 

The outline of this paper is as follows. In Sec. ||, we construct the pseudofermion action for 
one-flavor Wilson fermion, as well as that with the Hasenbusch acceleration [Q]. In Sec.[3|, we 
construct the pseudofermion action for one-flavor domain-wall fermion. In Sec.Q, we construct 
the pseudofermion action for one-flavor optimal domain-wall fermion. In Sec. ^ we compare the 
HMC simulations of 2 degenerate flavors with that of one flavor, for lattice QCD with optimal 
domain-wall fermion. In Sec. ^, we conclude with some remarks. 



2. Wilson fermion 

In this section, we derive a positive-definite pseudofermion action which gives exactly the 
determinant of the Wilson-Dirac operator 

D w = W + m + Ify » = ( <* + " > ' £ * °» V (2. 1, 



where 

~ 2 



W = [U^S^y + ulix- il)8 x ^ y ] +4, tp = \ [f/ /i (x)5 JC+Aj ,-f/^(x-M)4-A,y] . 



and Oi (i = 1,2,3) are Pauli mattices. Using the Schur decomposition, the determinant of (2.1) can 
be written as 

detD w =det(W + m) 2 detW//, (2.2) 
where W~h is the Schur complement of Dw, i.e., 

W H = W + m-Yt tl —^—t v ala v . (2.3) 
f^ v W + m » 

The positive definiteness of Wh is asserted as follows. For any background gauge field, the 
eigenvalues of W and (a -t) satisfy the inequalities: < X(W) < 8, and |A(^„(T^)| < 4. It 
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follows that Wh is positive-definite for m > 4. Now, if m is decreased from 4 to a smaller value, 
then the smallest eigenvalue of (W + m) is also decreased. Suppose that at m = m cr (the largest 
value), it becomes zero, then det(W + m cr ) = 0, and Wff(m cr ) becomes singular, thus Wh is no 
longer positive-definite. From ( |2.2| ), we see that whenever W# (ra) has a zero eigenvalue, Dw(m) 
also has a zero eigenvalue. On the other hand, when Wh (m) is singular due to a zero eigenvalue of 
(W + m), Dw(m) does not necessarily have a zero eigenvalue, which follows from the inequality, 
h m j n (W + m) < Re(X(Dw(m))) < X max (W +m). In other words, the largest m (i.e., m cr ) which 
gives det(W + m) = is equal to or larger than the largest m which gives detDw(m) = 0. 
Thus the pseudofermion action for one-flavor Wilson fermion can be written as 



S PF =<&](W + m)- z <& l +®l(W H y i ®2=®](W + m)- z <& l -(0 $>l)H w v II (2.4) 

where Hw = JsDw, <I>i is a pseudofermion field without Dirac index, and <I>2 is a pseudofermion 
field with 2 spinor components. 

To generate the pseudofermion field <I>2 from a Gaussian random noise field S2, we need 
to take the square root of Wh, i.e., <I>2 = ^Wh^-2 ■ Here we use the Zolotarev optimal rational 
approximation for the square root, 

O2 = V^S 2 ^ (pq + £ -P- J — ] S 2 , (2.5) 

\ i=\qi + W + m-Z iiy t^ WT ^t v ol 1 o v J 

where po, pi and q\ are expressed in terms of Jacobian elliptic functions. At first sight, the opera- 



tions in (2.5) look formidable. However, since Wh is the Schur complement of Dw, each term in 



(2.5) can be obtained by the inversion of (Dw +qiP~), i.e., 




P- [PO + £ Pl (Dw + q,P-y l L • (2-6) 



Note that one cannot apply the multi-shift conjugate gradient algorithm in ( |2.6P , since P_ does not 
commute with D^- However, these A^ app number of inversions can be speeded up using chronolog- 
ical inversion method [S]. 

2.1 Wilson fermion with Hasenbusch acceleration 

The idea of Hasenbusch's method [Q] is to introduce a heavy pseudofermion field (with mass 
M >m) such that the speed of simulating both det(Diy(ra) /Dw(M)) and det(Dw(M)) is faster than 
that of det(D w (m)). To apply Hasenbusch's method, we consider 

detfM^ =detf^VdetfBM^ (2 . 7) 



where 



W H (m)=W + m-Tt^— — Mfy<4- ( 2 -8) 
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Now, using the Schur decomposition, we can prove that 

det (D w {0)+miP+ + m 2 P- ) = det ( W + m i ) 2 det [W H (m 1 ) + A m ] = det (W + m 2 ) 2 det |F H (m 2 ) - A m ] , 

(2.9) 

where A m =m2-mi > 0. Using this relation, we can rewrite ( |2.7[ ) as 



det f ^pl) = det f ^ ) det f^M~M 
\DwMJ VW ff (m 1 )+A w 7 V W H {m 2 ) J 



Next we assert the positive definiteness of both operators on the RHS of ( |2.1C ). First, we assume 



that m.2 is sufficiently large such that (W -\-ni2), Wjiimi), and W uijni) are positive-definite. Then, 
the second factor on the RHS of ( |2.10D is positive-definite for m\ > m\ where mf is the largest 
value satisfying det[D^(0) + mfP + + /M2P-] = 0, as a consequence of ( |2.9| ). Now, for the operator 
(mi)(Wff (mi) +A m )~ 1 = (1 + A OT /Wff(mi)) , we note that even if W# (mi) is not well-defined, 
W H (mi)- 1 =P-H w ( mi )- l P- can still be well-defined. From (|2TT0|) , det ( W H (m 1 ) ( W tf (m 1 ) + A m ) ~ 1 ) 




is positive-definite if both detDvp(mi) and det[W#(m2) — A m ] are positive-definite. It follows that 
both operators on the RHS of ( |2.1CJ| ) are positive-definite for mi > max (mf ,mf), where mf has 
defined above, and rrv\ is the largest value satisfying det Dvp(mf) = 0. 

Thus the pseudofermion action for one-flavor Wilson fermion with Hasenbusch acceleration 
can be written as 

^ W H {mi J 

*;* 3 + *;* 4 - A m (o *y h w i (mi ; 

where *3 and *4 are pseudofermion fields with 2 spinor components. 

3. Domain-wall fermion 

The basic idea of domain-wall fermion is to use N s layers of 4-d Wilson-Dirac fermions (with 
nearest neighbor coupling in the 5-th dimension) such that the exactly chiral fermion fields emerge 
at the boundary layers in the limit N s —> °° [^]. The domain- wall fermion operator can be defined 

as 

#dwf(m) = W-m + £ y M f M +M(m) (3.1) 
M(m) = P + M + {m)+P-M-{m) (3.2) 

„ , s f^' )S -5^_i, 1<j<A^, [ 5,'..v-5 $ ', 4+ i, 1<j<A^, 

Af+(m)^/ = ^ M_(m) s/ = (3.3) 

[ cV,. v + mSy^ ,j = l, [ 5,',.v + m5 s ',i , 5 = N s , 

where m is the fermion mass, and tuq £ (0,2) is a parameter called "domain-wall height". Using 
the Schur decomposition, the determinant of domain-wall fermion operator can be written as 

det^ dwf (m) =det[W-m + M + (m)] 2 det#//(m) = det [IV - m + M_(m)] 2 det W H (m), (3.4) 
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where 



W H {m) =R 5 \ W-m +M^(m)-t^ 
W H (m) =R 5 [W -m + M+(m)-t^ 



1 



tvolo 



W -m +M + ( m y v * 
1 



-tv<5n<5 v 



(3.5) 



(3.6) 



"W — m + M-(m) 

Here R$ is the reflection operator in the 5-th dimension, (Rs) ss i = S s jy s+S r^i, which is introduced 
such that Wn{m) and Wijitn) are hermitian. After incorporating the contribution of the Pauli- 
Villars fields, the fermion determinant of domain-wall fermion becomes 

det ^ dwf (m) det [W - m + M+ (m) ] 2 det % (m) 



det dwf ( 1 ) det [W - mo + M+ ( 1 )] 2 det W H ( 1 ) 
Using the Schur decomposition of [^dwf(l) — (-^+(1) ~M+(m))P+], we obtain the relation 



det[W-m + M + (m)] 2 -det[#ff(m)+A_(m)] = det [W - m Q + M_(1)] 2 • det 



where 



7 tf(l)-A H 



[ A +H]^' = (1 -m)8s,N s 8 s i, Ns , [A_(m)] s/ = (l-m)S sA 8^. 
Using (3.8), we can rewrite ( |3~77| ) as 

det [W H {m)/{W H (m) + A_ (m) )] • det [(W H {\) - A + {m)) /W„(\)] 



(3.7) 

(3.8) 
(3.9) 



and its inverse 



det 



1+A_(m) 



%(m) 



•det 



l+A+(m) 



1 



/T H (l)-A + (m) 



can be used to construct the pseudofermion action. Using (f3.9|), we can simplify (3.1C) to 



det 1 + (1 



m) 



1 



W H {m) 



s=s'=\ 



•det 1 + (1 



1 



[m 



(3.10) 



(3.11) 



s=s'=N, 



Thus we can write the pseudofermion action for one-flavor domain-wall fermion as 

S PF = - (1 -m) (0 *J) te/? 5 ^dwf(m)]; = 1 y = i 





+ ^<I> 2 + (1 - m) (*| 0) [y 5 R 5 ^dwf(l) - A+(m)P + ]; = V = 



where and <J>2 are pseudofermion fields (on the 4-dimensional lattice) with 2 spinor components. 
Now we assert that the operators in ( J3- 11 ) are positive-definite for m > 0. At m = 1, they are equal 
to the identity operator, thus are positive-definite. As m is decreased to be less than one, their 
positive definiteness will lose only if either of the determinants in ( 3.11 ) becomes zero or singular. 
Using [W,M_(1)] = \W,M + (m)] = 0, and the fact that the eigenvalues of M_(l) and Af+(m) have 
non-zero imaginary parts, we immediately see that (W — niQ + M_(1)) and (W — mo+M + (m)) 
cannot have zero eigenvalue for m > 0. Thus the operators in (3.11) are well-defined for m > 0. 
Furthermore, since (3.10) is equal to the determinant of the massive overlap-Dirac operator with 
polar approximation for the sign function of H = y$D w (2 + D w ) , thus for m > 0, it must be 
positive and it follows that the operators in (3.1 1) are positive-definite. 
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4. Optimal domain-wall fermion 

For finite N s , the chiral symmetry of the domain-wall fermion is not exact. However, one can 
attain the (mathematically) maximal chiral symmetry by assigning each layer in the 5-th dimension 
a different weight (O s , according to the formula derived in Ref. ^7§. Then the effective 4D Dirac op- 
erator is exactly equal to the overlap-Dirac operator with Zolotarev optimal rational approximation 
to the sign function of H w . The optimal domain-wall fermion operator ^ is defined as 

%p,{m) = (QD w (-m )(l+L(m)) + (l-L(m)), L(m) = P + L + (m) + P_L_(m), (4.1) 

where 

T , v fe-i, Ks<N„ jS s > tS+h \<s<N s , 

[ -mb s i^ s , s = 1, [ -md s ^i , s=N s . 

Here ft) denotes the diagonal matrix of the weights {(O s , s = 1, • • • ,N S }. Moreover, we can construct 
ft) such that it satisfies R5C0 = (0R5. To derive a positive-definite pseudofermion action for one 
flavor, we can start from 

@' op t(m) = ft)D w (-mo) + (1 -L(m))/(1 +L(m)) (4.3) 

which is different from 3$ op t{m) by a matrix factor independent of the gauge field. Using the 
same technique we have used for the case of domain-wall fermion, we obtain the positive-definite 
pseudofermion action as follows 1 



S PF = +<I>^ 2 - ([^) (0 (-1)'*}) [ Y5 R 5 9' opt (m)]-l, 

1-1.(1) 1-L(m 



+ (its ) «-"^ o) -r, - ^ )p + V T'\ 

where and <I>2 are pseudofermion fields with 2 spinor components on the 4-dimensional lattice. 




5. Numerical test 

We compare the efficiency of HMC simulation of 2-flavor and ( 1 + 1 ) -flavor QCD with optimal 
domain- wall quarks, on the 12 3 x 24 x 16(N S ) lattice. For the gluon part, we use Iwasaki gauge 
action at )3 = 2.30. In the molecular dynamics, we use the Omelyan integrator and the Sexton- 
Weingarten multiple-time scale method |^]. The time step for the gauge field (AToauge) is the same 
for both 2-flavor and (1 + l)-flavor cases, while the time step (Atpp) for the pseudofermion field in 
the (1 + 1) -flavor case is 4 times larger than that for the 2-flavor case such that the acceptance rate 
is roughly the same for both cases. We use conjugate gradient (CG) with mixed precision for the 
inversion of the quark matrix (with even-odd preconditioning). The length of each trajectory is set 
to 2. After discarding 300 trajectories for thermalization, we accumulate 100 trajectories for the 
comparison of efficiency. Our results are given in Table.[j]. We see that acceptance rate is almost 
the same for (1 + l)-flavor and 2-flavor simulations. If auto-correlation time is the same, then the 
efficiency of HMC can be estimated by the total acceptance divided by the CG iteration number, 
and the efficiency ratio for 2-flavor and (1 + l)-flavor is about 3 : 2. 

1 Detailed derivation will be given in a forthcoming paper. 
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m 


N t 


< Bj /io 3 


< u, /io 3 




Acceptance 


Acceptance /N^^ 


0.019 


1 + 1 


75(1) 


260(1) 


345(1) 


0.88(3) 


2.6(1) x 10~ 6 


2 


0.6(1) 


239(2) 


240(2) 


0.90(3) 


3.8(2) x 10~ 6 


0.038 


1 + 1 


54(1) 


125(1) 


179(1) 


0.90(3) 


5.0(2) x 10" 6 


2 


0.6(1) 


112(1) 


113(1) 


0.91 (3) 


8.0(3) x 10" 6 



Table 1: Comparison of HMC efficiency for the 2-flavor and (1 + l)-flavor QCD with optimal domain-wall 
quarks. The step size for the gauge field ATcauge is 0.007(0.010) for m = 0.019(0.038). while the step size 
ATpp for (1 + l)-flavorpseudofermions is 0.14(0.20) for m = 0.019(0.038), which is 4 times larger than that 
for the 2-flavor case. Here, n2^\ N^f \ and are the average CG iterations for one trajectory (for 

generating initial pseudofermion fields, molecular dynamics, and their sum respectively). 



6. Concluding remark 

We have constructed a positive-definite pseudofermion action for one fermion flavor in lattice 
gauge theory, with determinant exactly equal to that of the lattice Dirac operator, for the Wilson- 
Dirac operator, and the (optimal) domain-wall fermion operator respectively. Evidently, these one- 
flavor positive-definite actions are useful for the HMC simulations in nonperturbative studies of the 
Standard Model (SM), as well as those beyond the SM. Our numerical tests show that the step size 
Atpp (for the pseudofermion field in the molecular dynamics) for the (1 + l)-flavor QCD can be 4 
times larger than that for the 2-flavor case, yet with the same acceptance rate. The efficiency ratio 
for HMC with 2-flavor and (1 + l)-flavor is about 3 : 2. With the positive-definite pseudofermion 



action (4.4) for the strange quark, TWQCD Collaboration is performing the HMC simulation of 



(2+l)-fiavor QCD on the 16 3 x 32 x 16 lattice, using a GPU cluster [[Kj]. 
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